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Modelling of energy flow in forest decomposer communities 
(Zur Funktion der Fauna in einem Mullbuchenwald 19) 


David H. Wise 


Synopsis 


Ein vorläufiges Modell für den Energiefluß in der Zersetzer-Gemeinschaft eines Waldbo- 
dens wird durch Simulationen auf Schlüssigkeit geprüft. Das Modell besteht aus drei Steu- 
ergrößen (‘driving variables’; Laubstreufall der Buche, Laubstreufall anderer Baumarten, 
Temperatur der Streuschicht und des Bodens) sowie neun dynamischen Zustandsvariab- 
len (‘dynamic state variables’; Laubstreu, Saprophage und Zoophage). Raten der Verän- 
derung der Zustandsvariablen werden durch nichtlineare Differentialgleichungen beschrie- 
ben. Werte für maximale Konsumptionsraten, Assimilationseffizienz, Respirationsraten 
und Korrekturfaktoren für Temperaturabhängigkeit wurden der Literatur entnommen. Als 
Lösung in seiner ursprünglichen Konfiguration ergab das Modell für die Zustandsvariablen 
Vorhersagen der Dichte, die nahe an den Durchschnittswerten für einen Kalkbuchenwald 
bei Göttingen lagen. Das Modell war jedoch sehr empfindlich gegenüber geringfügigen 
Änderungen in Parametern des Energieflusses, eine Eigenschaft, die die Anwendung des 
Modells in seiner jetzigen Form einschränkt. Diese Sensitivität wird diskutiert im Hinblick 
auf den verwendeten Modellansatz und auf die Entwicklung eines realistischeren Modells 
für den Energiefluß in der Zersetzerkette von Wäldern. 


Computer simulation, energy flow, forest decomposer communities, mathematical model- 
ling, sensitivity analysis. 


1. Introduction 


Ecologists increasingly have advocated mathematical modelling as a tool to summarize, 
examine and develop hypotheses about complex ecological systems. WIEGERT (1979), 
for example, has proposed that mathematical modelling of population processes can im- 
prove our understanding of ecosystem phenomena such as energy flow and nutrient cyc- 
ling. Two major types of systems models are used in ecological research: 


- Multiple regression models. Their primary goal is prediction of a particular system's futu- 
re state from empirically confirmed past patterns. Regression models do not necessarily 
incorporate explicit hypotheses about the processes underlying the correlations be- 
tween variables that comprise the model. 


- Functional models. The motivation behind their development is uncovering unforeseen 
implications, or consequences, of carefully specified ecological hypotheses. Under- 
standing the relationships between basic processes, rather than generating accurate 
site-specific predictions, is the primary objective. As with regression models, comparing 
the behaviour of real systems with the model’s predictions constitutes the main avenue 
of validation. 
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Functional modelling is the focus of the research addressed in this paper. This approach 
to modelling complex ecological systems is relatively recent and has not been done fre- 
quently with litter-soil communities. However, a few ecologists have attempted to model 
such decomposer systems. GIST and CROSSLEY (1975) constructed a linear, ten-com- 
partment model to describe the flow of calcium and potassium through the invertebrate 
food web of a leaf-litter community. Radiotracer studies provided estimates of rates of ex- 
change between some compartments. Data on feeding habits and the influence of tempe- 
rature on activity and reproduction were taken from the scientific literature. Simulation stu- 
dies with the model revealed that, in terms of their influence on rates of flux of these two 
nutrients, collembolans and cryptostigmatid mites were the most important saprovores, 
and spiders and mesostigmatid mites were the major predators in this litter community. 


Gist and Crossley’s model consists of linear differential equations - the simplest ones to 
handle mathematically. O'NEILL (1971) and WIEGERT (1975) argue that non-linear equa- 
tions more realistically represent the biological processes in decomposer systems. An ex- 
ample of this approach is the sub-model describing the decomposers in the terrestrial eco- 
logical energy model (TEEM) of the eastern deciduous forest biome section of the US IBP 
(SHUGART et al. 1974). This model describes energy transfers within a Liriodendron fo- 
rest. In the decomposition sub-model the fauna account directly for only about 5% of the 
loss of energy from the forest floor. However, computer simulations with the model sug- 
gest that removing the fauna from the system would produce a 20% increase in the 
amount of standing litter. Generation of such non-intuitive predictions is one of mathemati- 
cal modelling’s strengths as a research tool. Shugart and his co-authors point out that be- 
havior of the TEEM decomposition sub-model agrees with the outcome of a field manipu- 
lation by WITKAMP and CROSSLEY (1966), in which destruction of much of the fauna 
with naphthalene produced a 25% increase in the amount of litter. 


KOWAL (1971) modelled the interactions within the food web on the floor of a pine forest 
with a set of non-linear differential equations that incorporated reciprocal effects of preda- 
tors and prey upon changes in population size. The driving (input) variables were air tem- 
perature, net precipitation and the rate of input of pine litter. The dynamic state variables 
i.e., the variables that define the state of the system and whose rates of change with re- 
spect to time are described by differential equations, were litter moisture, amount of litter 
and biomass of nine faunal groups: spiders, ants, centipedes, three categories of mites 
and three categories of collembolans. Kowal made a preliminary examination of some of 
the model’s properties, but he did not compare its solution with the behavior of an actual 
decomposer food web. 


| have used Kowal’s basic approach as the starting point for my modelling effort, but have substanti- 
ally modified the non-linear differential equations in order to incorporate more detailed aspects of pre- 
dator-prey interactions. 


2. Overview of the model 


The model is a preliminary version that is being developed into amore complete model of 
energy flow through decomposer communities inhabiting beech forests on limestone, such 
as the extensively studied Göttinger Wald. Available data on faunal densities, rates of litter 
fall and soil/litter temperatures from this forest served as the focus for model development 
(SCHAEFER 1982, 1989 and unpublished data from several researchers). Nevertheless, 
the goal of the first phase of the modelling project is not to model this particular decompo- 
ser community. Rather, | wish to examine the properties of a simple model in order to de- 
termine whether a modelling approach that relies primarily upon simple hypotheses about 
energy flow can explain patterns of change in the abundance of major categories of or- 


328 


cianisms in such a system. The conclusions obtained so far relate to properties of the mo- 
del itself. The equations are not yet sufficiently refined to make inferences from the mo- 
del’s behavior about the importance of particular paths of energy flow in the Göttinger 
Wald or similar forests. 


Tab. 1: Model Variables. 


Driving (Input) Variables: 


D(1) --- Daily rate of input of beech leaves 

D(2) --- Daily rate of input of non-beech canopy leaves 
and herbaceous vegetation 

D(3) = Daily temperature of soil-litter interface 


State Variables: 


Primary Resources: 


X(1) Fresh beech litter and associated microflora 
X(2) Edible litter (non-beech and aged beech litter) 
and associated microflora 
Saprovores: 
X(7) Lumbricids 
X(9) Sciarid Diptera 
X(10) Diplopods and isopods 
X(11) Collembola 
Predators: 
X(21) Predaceous Coleoptera 
X(22) Chilopods 
X(24) Spiders 


This preliminary model consists of three driving variables and nine dynamic state variab- 
les (Tab. 1). The non-consecutive numbering of the variables highlights the models flexi- 
bility for future research, in that the FORTRAN program coded for the model allows up to 
26 variables to be incorporated. In this preliminary version some major groups have been 
omitted, primarily for logistical reasons associated with initial simulation studies. Four sa- 
provore categories consume one category of litter [X(2)] and three categories of predator 
consume either two, three or four saprovores (Fig. 1). The only exception to this pattern is 


LUMBRICIDS 
x(7) 


PREDACEOUS 
COLEOPTERA 


IN, 
CHILOPODS 
X(22) 


SPIDERS 
X(24) 


FRESH BEECH 
LITTER 
X(1) 


SCIARID 
DIPTERA 


EDIBLE 
LITTER 
X(2) 


DIPLOPODS & 
ISOPODS 


COLLEMBOLA 
X(11) 


Fig. 1: Pattern of energy flow in the model. 
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Tab. 2: Model parameters 


EXTCR 


WTHFR 


SW; 


ASij 

Rj 
TC,TR; 
TS; 


DIM, 


Extinction criterion: 
% of initial value of state variable below which the value is set = 0 
for the remainder of the simulation 


Fraction of beech litter converted to edible litter 
category each September 


Prey selectivity factor: 
relative intensity of feeding upon resource i by consumer j 


Maximum consumption rate (kJ/kJ/d) of consumer j 
Density of resources at which consumption rate = .5C; 


Predator-switching criterion: 
% of resource i in diet of consumer j below which j does not feed oni 


Assimilation efficiency of j feeding on resource I 
Respiration rate (kJ/kJ/d) of state variable j 
Temperatures at which Cj and Rj were measured 
Temperature sensitivity of Cj and R; 


Density-independent mortality of consumer j 


that in one of the simulations of the model’s behavior, a small amount of energy flows di- 
rectly from fresh beech litter [X(1)] to the lumbricids [X(7)] (indicated by a fainter line con- 
necting these two state variables in Fig. 1). 


The parameters of the model appear in Tab. 2. Maximum consumption rate (Cj) and Hj, 
the density at which the actual consumption rate equals .5C;, define the functional respon- 
se (Fig. 2). The general equation describing the rate of change in consumer state variable 
j that utilizes resource i and is eaten by state variable k is: 


axat = D ATTY CCI - Je Ge re 


1 


= Ry(e >) (D(3) - TRJ) y xq) 


“KZ X(n), (TS, (D(S) - TC) 


) X(k) 


( 


i 


— DIM, 


The attack function, ATTij, determines the relative proportion of resource i in the diet as a function of 
the switching criterion, the prey selectivity factor, and the proportion of i in the potential resource 
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spectrum of consumer j. The differential equations describing the growth of the two primary resour- 
ces, fresh beech litter [X(1)] and edible litter [X(2)], differ from this general equation in the first term. 
For X(1) this term is D(1), and for X(2) it is {D(2) + WTHFR[X(1)]}. The remainder of both equations is 
the same as that for the consumers, except that for X(1) the additional amount WTHFR[X(1)] is sub- 
tracted each September. 


The initial value of each state variable was set equal to its average in the Göttinger Wald 
ior the period April 1981 - March 1982. Daily rates of input of beech litter, D(1), were 
based upon rates for the Göttinger Wald from June 1981 through December 1983. Chan- 
ges in the second driving variable, D(2), reflected the production of non-beech canopy lit- 
ter from October 1981 through December 1983, and the annual pattern of production of 
herbaceous litter in the Göttinger Wald. Temperature input, (D3), was based upon the 
daily average at the soil/litter interface from April 1981 through January 1986. Available 
data for the driving variables do not span periods of time long enough for useful solutions 
ofthe model. Therefore, the above data were used to calculate an average value for each 
day of the year; these average yearly patterns for the three driving variables were then 
used as repeated input in order to simulate the behavior of the model for up to 25 years. 


Each solution to the set of differential equations for a particular set of parameter values is 
an integral curve obtained by the Euler approximation with a time increment of .5 day. In 
the examples that follow, each solution for a particular state variable is plotted along with 
the average value of that variable in the Göttinger Wald (the horizontal dashed line in the 
figures). 


3. Development of the standard solution 


For the standard solution and all subsequent modifications of the model parameters, 
EXTCR = .01, WTHFR = 1, SWij = .01 and DIM; = 0. In the standard solution Sij = 1 for all 
energy pathways and it is assumed that lumbricids do not feed upon freshly fallen beech 
litter (i.e. S17 = 0). Saprovores and predators have different functional responses in the 
model, reflecting different hypotheses about the extent of food limitation of each. Values of 
Hj for saprovores were set so that the actual consumption rate of each = .99C; at a resour- 
ce density equal to .1 of the sum of the initial values of all the resources consumed by that 
saprovore state variable. Unlike saprovores, the predators in the model display a pronoun- 
ced Type II functional response (Fig. 2); Hj equals the sum of the initial values of all re- 
sources consumed by j. 


a 


kJ resource/ 
kJ consumer/ 
day 

(Y) 


Y=c(1-e*%) 


Where: 
k = -(1/H) In(.5) 


kJ resource/sq. m 


(X) 


Fig. 2: Form of the functional response in the model. 
(1) resource density does not affect feeding rate at most resource densities -- characteristic 
of saprovores in this model. 
(2) typical type Il functional response curve -- characteristic of predators in this model. 


331 


Variable Input Constant Input 


X(1)-FRESH BEECH LITTER 


8 0 2 4 6 8 
X(7) — LUMBRICIDS 


8 o 
X(9)- SCIARID DIPTERA 


YEARS 


Fig. 3: Model solutions, over 8.3-year interval, using standard parameters (Tab. 3) and 
two types of inputs. 
Variable input: variation in driving variables reflects daily variation in Gôttinger Wald. 
Constant input: constant values of driving variables, equal to yearly averages. 
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Variable Input | Constant Input 


X(10)- DIPLOPODS & ISOPODS 


0 
xt21)- PREDACEOUS COLEOPTERA 


Fig. 4: 
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X(7)- LUMBRICIDS X(9)- SCIARID DIPTERA 
300 
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X(10)- DIPLOPODS & ISOPODS X(11)- COLLEMBOLA 


0 6 12 18 24 0 6 12 18 24 
X(21)- PREDACEOUS COLEOPTERA X(22)- CHILOPODS 


6 
X(24)- SPIDERS 


YEARS 


The standard solution for saprovore and predator state variables. 

X(1) and X(2) are omitted from this and following figures because their patterns of fluctua- 
tion do not change substantially over longer time intervals nor with the parameter perturba- 
tions reported here. 


values of Cj, ASij, Rj, TCj, TRj and TS; were taken, whenever possible, from the published 
lierature (DOEKSEN, DRIFT 1963; ANDERSON 1970; KOWAL, CROSSLEY 1971; 
McBRAYER, REICHLE 1971; RUNGE 1973; HEAL, MACLEAN 1975; PERSSON, LOHM 
1977; REICHLE 1977; HUMPHREYS 1979; ALBERT 1983; SATCHELL 1983). In many 
cases several values were available, but for some parameters an educated guess was ne- 
cossary. Some of the values originally selected were adjusted so that each saprovore 
“oubled in size after a year in the absence of predators and other saprovores, and with a 
superabundance of X(2). In addition to giving the predators a pronounced functional re- 
sponse to changes in resource abundance, food limitation was modelled further by adju- 
sing Rj of the predators so that the actual rate of consumption = RjASi; when resource 
densities equalled the average values for the Göttinger Wald. 


Respiration rates of the litter/microflora complex were chosen so that in model simulations 
with no saprovores, the rate of disappearance of X(1) and X(2) mimicked the observed ra- 
te of disappearance of leaves from fine-mesh litter bags placed in the Göttinger Wald for a 
year (unpubl. data). Rı and R2 were the only parameters that were based upon the obser- 
ved behavior of state variables in the Göttinger Wald; all other parameter values came 
from the literature or reflected particular hypotheses about interactions within decomposer 
communities. The values of the basic parameters of energy flow for all the state variables 
are listed in Tab. 3. 


This combination of parameters yielded a solution of the model over an interval of 8 years 
‘hat agreed satisfactorily with the average value of each state variable in the Göttinger 
Wald. This solution resulted from input that reflected daily variation of the driving variables 
in the Göttinger Wald (actual values during the first few years, followed by the average ye- 
arly pattern repeated for the remaining number of years). An initial goal, however, was to 
study the model’s behavior in the absence of environmental fluctuations, in order to reveal 
its structural properties. In particular, | wanted to investigate the consequences of the 
strong species interactions that are a key point of this model. Therefore, the model was 
solved over the same time interval using constant input; each driving variable equalled its 
yearly average [D(1) = 17.3, D(2) = 2.5, D(3) = 6.0]. The solution with constant driving va- 
riables is remarkably similar to that based upon the variable input that reflects natural en- 
vironmental fluctuations (Fig. 3). 


Because of this close agreement, subsequent simulations were done with constant driving 
variables and the solution interval was expanded to 25 years. The behavior of the model 
under these conditions, with parameter values equal to those described above, is termed 
the standard solution (Fig. 4). The standard solution was judged to mimic adequately the 


Tab. 3: Values of parameters in the standard solution. 
Parameter symbols are explained in Tab. 2. Absent from this table are parameter values gi- 
ven in the text: SW, S and DIM. 


Para- State Variable 
meter X(1)  X(2) AN) X9)  X(10) X01)  X(21)  X(22)  X(24) 


C ae = 11 .06 .04 .09 .03 .05 .05 
H --- = 50 50 50 50 212 223 22 

TC --- es 10 10 10 10 10 10 10 

R .00014 .00379 .0025 .0015 .005 .022 .012 .02 .015 
TR 6 6 10 10 10 10 10 10 10 

TS .069 .069 .09 .069 .10 .09 .09 11 .14 
AS -- --- .05 .07 co 36 .80 .80 .60 
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behavior of a real decomposer community because the integral curve of each state vari- 
able is close to its average value in the Göttinger Wald: differences between predicted 
and average values are well within variations in empirical values due to actual fluctuations 
and sampling error. The marked annual oscillations in densities of the saprovores reflect 
seasonal variation in the availability of edible leaf litter, X(2), due to changes in litterfall 
D(2) and the movement of non-edible fresh beech litter [X(1)] into the edible-litter category 
each September. The oscillations with longer periods likely reflect the destabilizing effects 
of the predator-prey interactions in the model. 


4. Some preliminary simulation experiments 


A major goal of constructing functional mathematical models is to examine consequences 
of the hypotheses that compose the model. Results of research on the model can point to 
potentially fruitful areas of empirical research. Two examples of such simulation experi- 
ments are described next. 


4.1 Feeding of lumbricids on fresh beech litter 


Lumbricids most likely do not restrict their feeding to non-beech litter and beech litter that 
has aged a year or more, as occurs in the standard solution. What are the consequences 
for the entire system if lumbricids feed on fresh beech litter, X(1), in addition to X(2)? The 
answer, based upon the behavior of the model, is quite surprising. Allowing lumbricids to 
feed upon freshly fallen beech leaves at 1% of the intensity with which they feed upon the 
preferred resource (S17 = .01, S27 = 1) causes the diplopods and isopods [X(10)] to grow 


X(7)- LUMBRICIDS X(10)— DIPLOPODS & ISOPODS 


2 


m 


0 6 12 18 24 0 6 12 18 24 
X(21)- PREDACEOUS COLEOPTERA X(22)- CHILOPODS 


kJ/ 


YEARS 


Fig.5: Solution of model with S17 = .10 (S17 = 0 in standard solution). 
Representative patterns for four state variables are shown. Compare with standard solu- 
tion in Fig. 4. 
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exponentially and drives the spiders to extinction within 7 years. Even more dramatic ef- 
tects result from increasing Sı7 to 10%. The spiders are still the only organisms to go ex- 
tinct, but the growth patterns of the other state variables change to pronounced cycles 
with amplitudes and periods quite different from any appearing in the standard solution 


(Fig. 5). 


These dramatic changes in the behavior of the model are non-obvious, non-intuitive con- 
sequences of altering a simple hypothesis about feeding relationships. The outcome of 
this simple simulation demonstrates the potential power of modelling as a tool for investi- 
gating hypotheses about complex ecological systems. If the model accurately represents 
natural processes, these simulation results would suggest that it is critical to our under- 
standing of the functioning of the entire community to know the fraction of freshly fallen 
beech litter in the lumbricid diet. 


4.2 Removal of predators 


Predators among the macrofauna might affect primary decomposition indirectly if they si- 
gnificantly limit densities of their saprovore prey. Whether saprovore densities are signifi- 
cantly depressed by predators is unknown, primarily because predator-removal experi- 
ments in natural decomposer communities are difficult to perform. One preliminary appro- 
ach to answering this question is to examine the behavior of models of decomposer 
communities from which predators have been removed. Excluding all three predator state 
variables from the standard model suggests that X(7) and X(9) are predator-limited, and 
that predation ameliorates the competitive effects of these saprovores upon X(10) and 
X(11) (Fig. 6). 


X(7)- LUMBRICIDS X(9)-SCIARID DIPTERA 


0 6 12 18 24 0 6 
X(10)- DIPLOPODS & ISOPODS X(11)- COLLEMBOLA 


kJ/m 


YEARS 


Fig. 6: Solution of model for the saprovores when all the predators have been excluded. 
Compare with standard solution in Fig. 4. 
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Because this model is a preliminary version, these two simulation experiments are pre- 
sented as examples, not as conclusions or predictions about processes occurring within 
any natural decomposer community. Before planning empirical studies to test further a 
models validity, one must scrutinize the model to determine if unreasonable assumptions 
lie hidden within its structure. A first approach is to examine the model’s sensitivity to 
some basic parameters, which for this model are those directly associated with energy 
flow i.e. assimilation efficiency, maximum consumption rate and respiration rate. 


5. Sensitivity of the model to changes in some basic parameters 


Results of sensitivity analyses of these three parameters are summarized below. In each 
perturbation the parameters for all consumer state variables [all state variables except 
X(1) and X(2)] have been displaced the same percentage from the value in the standard 
solution. Only one parameter is changed at a time i.e. values are set equal to those that 
generated the standard solution before the next perturbation is performed. 


A) 


MIDPOINT OF RANGE or "EQUILIBRIUM" DENSITY 


B) 


MIDPOINT OF RANGE or "EQUILIBRIUM" DENSITY 


Ê AS by 10% 7 = 0 RTE Pc by 10% 
20% 0 EI E ee 20% 
30% 0 Ele E | ¢ fle 100% 
| AS by 10% fioleiWbi le | © by 10% 
20% tt | o EJH Y| E 
30% t | E| E E | E 


MIDPOINT OF RANGE or "EQUILIBRIUM" DENSITY 


Respiration 
rate (R): 


$ R by20% 


0 : .65x -- 1.5x t : 1.5x = 
E : EXTINCT tt >5x H: <2x 


50% 
100% 
| R by 20% 


50% 


Fig. 7: Sensitivity of the model solution to changes in some basic parameters of energy 
flow. 
(A) assimilation efficiency. 
(B) maximum consumption rate. 
(C) respiration rate. 


Increasing the ASij by only 10% causes the average density of the Diptera to decrease by 
90%, the average density of the Collembola to increase over 8-fold, and the diplopods and 
isopods to go extinct. Decreasing ASij by 10% causes the Collembola and the spiders to 
go extinct. Changing AS; by only 30% causes even more dramatic changes in the aver- 
age densities of the state variables (Fig. 7A). 


Changing Ci produces patterns similar to those when ASj is perturbed (Fig. 7B). 
Somewhat surprisingly, the model appears to be less sensitive to changes in Ci; than ASij. 
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‚ncreasing Rij by only 20% causes the average value of one saprovore to increase over 8- 
told, the mean density of two predators to decline by 50%, and two state variables to go 
extinct. Decreasing Rij by only 20% produces four extinctions. Altering Rij clearly modifies 
the model’s behavior, and in addition, the patterns of the often dramatic changes are not 
always predictable (Fig. 7C). 


6. Conclusions and discussion 


Clearly the reasonably good behavior of the standard solution is misleading, because the 
model is excessively sensitive to minor changes in several basic parameters of energy 
flow. It is unlikely that the model as now formulated accurately represents the control of 
energy flow through a forest decomposer community. If the representation were valid, we 
would have to postulate that actual respiration rates, assimilation efficiencies, and maxi- 
mum consumption rates are so finely tuned, so carefully in balance, that any small shift 
would result in collapse of the system. To make such an assumption would be unreasona- 
ble. In nature rates change during growth and development of the organism, differ be- 
tween species within a state-variable category, and change in response to variation in 
physical variables. Furthermore, we cannot measure average rates with the degree of pre- 
cision that would be required to set correct parameter values for an empirical test of this 
model in a natural community. 


Such instability in a model is unreasonable and severely limits its utility. The extreme sen- 
sitivity to changes in respiration rates, assimilation efficiencies and consumption rates se- 
riously calls into question the advisability of modelling energy flow in decomposer com- 
munities with an approach that relies primarily upon using the energy-budget equation of 
an individual organism to describe the dynamics of pooled species populations. Even 
though differences in developmental stages have been ignored and species populations 
have been pooled, it is still possible that the basic approach utilized so far may be overly 
reductionist. It might be more fruitful to reduce the number of parameters that affect popu- 
lation growth and write simpler differential equations, adjusting the relatively few constants 
to agree with the results of field perturbation experiments that could be conducted with 
components of the decomposer community. Unfortunately, such field experiments are not 
always easily performed, and very few have been done to date. 


Despite the misgivings engendered by the sensitivity of the model’s parameters, it would 
be premature to reject completely the basic approach to modelling energy flow that has 
been used up to now. Possibly the sensitivity of the model can be reduced substantially by 
modifying its structure in one or more ways: (I) by incorporating a soil organic matter com- 
partment, a distinguishing feature of the food web of decomposer communities that may 
greatly affect community dynamics; (Il) by changing the pattern of weathering of beech lit- 
ter; (III) by adding state variables, both consumers and a greater variety of primary resour- 
ces; (IV) by incorporating additional driving variables, such as moisture content of the litter 
and soil; (V) by altering the pattern of feeding linkages; (VI) by modifying the resource 
switching criteria and the functional responses exhibited by all consumers; (VII) by incor- 
porating density-independent mortality. 


Mathematical models are, at best, caricatures of nature. They are valuable if they uncover 
counter-intuitive, non-obvious consequences of hypotheses, and if they aid in planning 
empirical studies. From this perspective mathematical models are never the end products 
of research, but are potentially powerful tools. So far their potential often appears to be 
more promise than reality, but certainly it is worthwhile to continue investigating their use- 
fulness for guiding research on complex systems. 
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